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O ' Abstract 

. . 

Q^' A simple and well-motivated explanation for the origin of dark matter is that 



it consists of thermal relic particles that get their mass entirely through elec- 



D . troweak symmetry breaking. The simplest models implementing this possibility 

predict a dark matter candidate that consists of a mixture of two Dirac neu- 
trinos with opposite isospin, and so has suppressed coupling to the Z. These 
^ ' models predict dark matter masses of m/jjvf ~ 45 GeV or m dm ~ 90 — 95 GeV 

' and WIMP-neutron spin-independent cross sections awiMP-n ~ 

pb. Current direct dark matter searches are probing a portion of the parameter 
space of these models while future experiments sensitive to crwiMP-n ~ 10~^ 
pb will probe the remainder. An enhancement of the galactic halo gamma ray 
and positron flux coming from annihilations of these particles is also expected 
across the ~ 1 — 100 GeV range. The framework further suggests an environ- 
mental explanation of the hierarchy between the weak and Planck scales and 
of the small value of the cosmological constant relative to the weak scale. 
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1 Introduction 

Stable particles produced thermally in the early universe with weak interactions and 
mass near ~ 500 GeV can account for the measured density of dark matter. This 
suggests a possible relationship between dark matter and electroweak physics. Conse- 
quently, most efforts to explain the origin of dark matter have focused on the Lightest 
Super-Partner (LSP) in the Minimal Supersymmetric Standard Model (MSSM) p]. 
This possibility is attractive because the presence of weak-scale supersymmetry may 
resolve the gauge hierarchy problem, predicts gauge coupling unification, and explains 
the connection between the dark matter scale and weak scale. Yet weak-scale super- 
symmetry was not discovered at LEPII nor have any conclusive indirect signatures 
been seen In light of this, the possibility that nature is seemingly fine-tuned 
has received renewed interest |2l EI , and so we reconsider a naive explanation for the 
existence of weak-scale dark matter that does not rely on supersymmetry: Can it be 
that the dark matter scale coincides with the electroweak scale because the dark matter 
particles get their mass entirely from electro-weak symmetry breaking? 

In this paper, we show that the answer to this question is yes, that such a frame- 
work is consistent with experimental data, and that it does not require elaborate 
fine-tunings. Moreover, it is surprisingly predictive and suggests resolutions to both 
the gauge and cosmological constant hierarchy problems within a landscape frame- 
work. 

In Section |21 we review the evidence for dark matter and briefly consider the 
experimental constraints on heavy (~ 20 — 100 GeV) neutrinos produced thermally 
as dark matter candidates. We conclude that, unlike neutrinos with invariant mass 
terms, a single neutrino with only a Higgs- Yukawa mass is always under- abundant. 
We note, however, that suppression of the neutrino coupling to the Z boson can 
remedy the many problems with a single heavy neutrino. 

In Section 121 we study the effects of Z suppression on the thermal heavy neutrino 
density. The suppression can be realized by adding to the Standard Model with a 
Higgs doublet a vector-like "dark sector" consisting of two lepton generations with 
opposite hypercharge and two neutral singlets ()3.1|) . A chiral symmetry prevents 
these states from mixing with standard model particles or acquiring invariant mass 
terms. After electroweak symmetry breaking, the stable ?7(l)£;M-neutral particle is a 
mixture of opposite-isospin neutral states, thereby suppressing its coupling to the Z. 
We present calculations of the dark matter abundance in this toy model ()3.2|) . 

In Section^ we revisit the model of 13.11 and another two- neutrino model allow- 
ing Z-coupling suppression, focusing on the radiative stability of the near-maximal 
mixing required to reproduce the measured ^Idm- The second model ()4.2|) is remi- 
niscent of the particle content of Split Supersymmetry |5j, with the gluino removed 
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and supersymmetric relations badly broken. We also point out in ()4.3|) that gauge 
coupling unification can be achieved at Mq ~ 10^^ GeV by extending the dark sector 
to include three of the vector lepton families of Sec. 13.11 

In Section El we consider the consistency of the toy models of Sections El and jH 
with current experimental constraints, and the prospects for detection in the near 
future in direct ()5.1|) and indirect ()5.2|) dark matter searches, as well as in colliders 
fl5.3|) . We focus on the region of parameter space that is most consistent with the 
observed dark matter density. We also consider the constraints placed on the model 
by precision electroweak observables (j5.4|) . 

Although our motivation is independent of the weak-Planck hierarchy problem, 
the models we consider do suggest possible resolutions to this problem within a land- 
scape framework. The crucial link in this argument is Weinberg's observation that 
structure formation requires that the cosmological constant be small jH] • By explicitly 
connecting the dark matter mass to the Higgs vev, the class of models suggested here 
relates the weak-A hierarchy to the weak-Planck hierarchy. We also further develop 
the suggestion in ^ that the gauge hierarchy may be explained using the structure 
principle because the existence of massive dark matter is tied to the stability of the 
Higgs vacuum. These two arguments are discussed in Section El 

2 Dark Matter: Evidence and Constraints 

The existence of dark matter is by now well established. The most direct evidence 
comes from comparing the luminous mass of a galaxy or cluster to its total mass, 
inferred by indirect means. Measurements of galactic rotation curves indicate that 
the mass of galaxies must be greater and more spread out than the luminous mat- 
ter they contain, suggesting a dark matter component Qdm > 0.1. Various means 
of estimating the mass of galaxy clusters, such as velocity dispersions and gravita- 
tional lensing, imply Qdm ~ 0.2 — 0.3 on the cluster scale. A less direct — but more 
precise — probe of dark matter is the anisotropy of the cosmic microwave background 
(CMB)j7]. Assuming a ACDM model (in which all unknown energy density is as- 
sumed to come from either a cosmological constant or non-relativistic dark matter), 
current measurements indicate a total matter component of f^M^^ = 0.134 ± 0.006, 
with a much smaller baryonic matter component of fi^/i^ = 0.023 ± 0.001 (which is 
consistent with the limits on Qb from the abundances of light elements produced in 
big-bang nucleosynthesis 0). Along with MACHO searches jHUHl! this is the principal 
observational evidence for the non-baryonic nature of at least some of the dark mat- 
ter. Finally, detailed modeling of our own galaxy suggest a local dark matter density 
of Pdm ~ 0.3 GeV/cm^. For more careful discussions of the evidence for dark matter 
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and its properties, we refer the reader to the reviews (TUl HH 112] • 

In hght of the increasingly precise evidence for a sizable dark matter component, 
many explanations of its nature and origin have been proposed (see |im ITU IT^ for a 
review of both particle physics-motivated and purely astrophysical candidates) . The 
elementary particle candidates are either (1) thermal relics, which begin in thermal 
equilibrium in the early universe but decouple from the plasma as the expansion rate 
of the universe exceeds their annihilation rates or (2) non-thermal relics, produced 
out of equilibrium (for example, from the decay of another particle that has itself 
decoupled from the plasma). The models we propose are of the first class. We further 
specialize to the case of dark matter particles that become non-relativistic before they 
decouple. 

In the hot early universe, the dark matter particles are kept in chemical and 
kinetic equilibrium with the surrounding particle bath by pair production, annihila- 
tion, and scattering. If their interactions with the bath are sufficiently strong, they 
remain coupled until after they become non-relativistic. As the universe cools below 
T ^ mciM, their density becomes exponentially Boltzmann-suppressed, quenching the 
annihilations. This "freeze-out" happens quite suddenly, typically at a temperature 
of Tf ^ . From this point on, the co-moving particle density remains approxi- 
mately constant. Thus, we can compare the observed dark matter density today to 
the predicted abundance, which is very roughly 



Assuming the maximum perturbative cross section ~ — for the dark matter 



particle's annihilations and requiring consistency of Equation ^ with the observed 
dark matter component places an upper bound on its mass, moM 300 TeV [12,. For 
a weak-interaction coupling, the mass scale of interest is near a TeV, coincidentally 
close to the scale of electroweak symmetry breaking. Weakly interacting particles 
with mass near a TeV (WIMPs) provide especially intriguing dark matter candidates 
because of their possible connection to electroweak physics. 

One simple WIMP candidate is a stable heavy Dirac or Majorana neutrino. Neu- 
trinos with a Dirac- or Majorana-type invariant mass have small enough annihilation 
rates to account for the observed dark matter abundance if they are heavier than 
~ 500 GeV, though this option seems to be ruled out by direct dark matter searches 
(see Sec. 15. Ij) . The best-motivated and currently well-studied WIMP candidate is the 
LSP in supersymmetric extensions of the Standard Model. In these models, the dark 
matter mass and electroweak symmetry breaking scale are close for a reason — because 
both are determined by the scale of supersymmetry breaking. 



0.1p6 X c 
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Another possible explanation of this proximity of scales is if, whatever may set 
the EWSB scale, the symmetry breaking itself sets the scale of the dark matter 
mass, which comes entirely from a Higgs Yukawa coupling. Neutrinos with Standard 
Model coupling to the Z that get mass from electroweak symmetry breaking are 
well ruled out by experiments. Z-pole width measurements at LEP have ruled out 
neutrinos lighter than and a mass above ^ would lead to over-annihilation. 
Thus QoAih"^ = 0.11 ± .011 cannot be obtained from a Standard Model-like neutrino 
of any mass. We revisit the possibility that the dark matter particle is given mass by 
electroweak symmetry breaking and observe that the problems associated with heavy 
neutrino dark matter can be resolved if the particle's coupling to the Z is suppressed. 



We consider a model suggested in ^ consisting of the minimal Standard Model, a 
lepton generation (Li, ^i) with the usual standard model charges, and one {L2,E2) 
with opposite hypercharge assignments. To this we add two SU{3) x SU{2)l x U{1)y 
singlets (si, S2) (this particle content is among the smallest anomaly-free extensions of 
the Standard Model in which all particles can obtain mass via EW symmetry break- 
ing). A chiral symmetry under which the Standard Model is neutral and the two 
new lepton doublets have charge opposite that of the four new singlets (or a discrete 
subgroup thereof) assures that all invariant mass terms vanish and mixing with Stan- 
dard Model leptons is eliminated, making the lightest particle exactly stable. This 
chiral symmetry plays a role analogous to that played by R-parity in supersymmetric 
models. 

The most general Lagrangian consistent with this symmetry takes the form 



where h"^ = —ia2h*, and implicit products of SU{2) doublets are taken to mean 
AB = A^i(j2B. Up to re-phasing, the 2x2 matrix K can be written as 



3 Dark Matter with Suppressed Z-Coupling 



3.1 A Toy Model 



Ant = Yih'LiEi + Y2hL2E2 + (/iLi, h'L2)K 




(2) 




(3) 



with V unitary. 

In terms of mass eigenstates, the Higgs- and Z-interaction Lagrangian for the 
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t^(l)£;M-neutral states takes the form 



h., ,x 9 ^ , ^Ti-^ „ ( cos(26') sin(26') \ ( s 

£ D 1 + - miVs + m'iV's' + ^ \ Z^{N, iV' . _ „ _ „ ^ ' 

V zcosb'iy \ sm(26') co?,[26) 



where 

'iV\ /A^i cos(^) - A/'a sm(0)\ f s\ 



(4) 

A^V V^i sin(^) + A^2 cos(^)y ' \s' ) "^^82}' ''^^ 
and ^ = ^-iid m' = -^k' are the neutral particle masses after h takes the vev 

{h) = -j=\ ) . The strength of the neutral fermions' coupling to the Z is suppressed 



y/2 \. 



by e = cos(2^) compared to that of a Standard Model neutrino. As we shall see, in 
the ^ ~ ^ regime, the particle is a viable dark matter candidate. Henceforth, we refer 
to the lighter neutral state as a Dirac fermion x, and the heavier as x'- 



3.2 Relic Density 

The framework for calculating relic abundances for thermally produced, stable par- 
ticles is discussed in Appendix |B| For a Standard Model Dirac neutrino, the relic 
density ignoring all co-annihilations is shown in Figure ^ together with the thermally 
averaged cross section {av) aX x = Xf. Reproducing the observed VLdm requires 
(cTf) ~ 2 pb c (double the estimate from [T] because both x aiid x contribute). As 
can be seen, Z-mediated annihilations into fermion final states dominate below the 
Z-pole so that the relic density coincides with the experimentally observed level only 
at ~ 5 GeV, a mass which has been ruled out by Z-width measurements at LEP. Even 
as the Z-mediated annihilation rate drops with increasing moM-, Higgs-mediated an- 
nihilations into ZZ and WW pairs are strong enough that the relic density never 
exceeds 10^^. Thus, a neutrino with Higgs- Yukawa-type mass can only be the dark 
matter if the Z-coupling is suppressed. 

The vector lepton toy model defined in 13. II allows this suppression as does another 
model introduced in Sec. 14.21 The Z-suppression in both cases results, as discussed 
above, from mass-mixing between gauge eigenstates with opposite isospin. There is 
always a corresponding enhancement of the Z-coupling between the heavy and light 
eigenstates, but if one is significantly heavier than the other this has little effect on 
the relic density. 

We have calculated the relic density of the lightest dark sector neutrino as a func- 
tion of its mass m^M and Z-coupling suppression factor e = cos 26*. Our calculations 
are based on the perturbative freezeout methods discussed in O CH CS] and the 
details of our calculation can be found in appendix ^ We ignore the effects of co- 
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Figure 1: Top: Relic density versus mass for a heavy Dirac Standard Model neutrino. 
Bottom: (av) at a; = 24 (near freeze-out) versus mass for a Dirac Standard Model 
neutrino. Both plots assume a Higgs mass of mh = 115 GeV. 
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annihilations, which are unimportant because the dark matter state can be made 
naturally much lighter than the other states that could co-annihilate with it. 

Figures El El andjUdisplay our results for Higgs masses of mh = 115, 140, and 160 
GeV respectively. In each case, we consider Z-suppressions e = 0.05, 0.1, 0.15, and 
0.2. The highest curves correspond to the lowest e, for which annihilation rates are 
lowest and the resulting relic densities highest. The viable neutrino mass below the 
Z pole is pushed up when the Z coupling is suppressed, from moM ~ 5 GeV for a 
Standard Model neutrino to mnM ~ 35 GeV with e < 0.1. As we discuss in 15.31 the 
Z- suppression of 0.1 and phase-space suppression from proximity of to mz/2 

is enough for a heavy neutrino in this mass range to have avoided detection at LEP, 
while further suppression is required to evade the null results of many direct WIMP 
searches. 

In the region m,£)M > there is a band of masses for which some e ~ 0.1 would 
produce the observed abundance. This band is bounded by the Z-pole and by the 
onset of annihilation into W^W~ and ZZ pairs. This region sits near the s-channel 
Higgs pole, so it is also chopped off where the efficient annihilation near the pole 
suppresses relic density. For a Higgs mass of mh = 115 GeV, the relic density ^hj^j^^ 
is consistent with experiment in the range of moM ~ 60 — 85 GeV. On the Z-pole and 
above this region, near ~ 90 GeV, the observed abundance can be reproduced 

by assuming very near-maximal mixing (e < 0.01). 

When Z-suppression is allowed, there is an interesting region of parameter space 
in which a heavy neutrino can reproduce the Qdm inferred from observations. The 
supression necessary is, however, far greater than one would expect from a generic 
choice of Yukawa couplings. In the following section, we examine the radiative sta- 
bility of this minimal mixing and its possible origins. 

4 Models With Maximally Mixed Neutrinos 

In this section, we discuss two minimal extensions of the Standard Model in which 
maximally mixed neutrinos with no invariant mass accounts for the dark matter. The 
first is the toy model discussed above. The second has a particle content reminiscent 
of Split-SUSY but with no supersymmetric relations among parameters. 

4.1 Vector Lepton Model 

As discussed above, the model of Section 13.11 generically mixes neutrino species of 
opposite isospin, resulting in supperssion of the Z-coupling. But near-maximal mixing 
(e < 0.15) is necessary to reproduce the observed relic density and to evade exclusion 
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Relic Density for mH=115, e=0.05, 0.1, 0.15, and 0.2 
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Figure 2: Top: Relic density versus dark matter mass for Z suppression of (top to 
bottom curves) 0.05, 0.1, 0.15, and 0.2. Bottom: Thermally averaged annihilation 
cross section {av) by final state at x = 24: (near freeze-out) versus mass, with Z 
suppression of 0.05. Both plots assume a Higgs mass of nih = 115 GeV. 
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Relic Density for mH=140, e=0.05, 0.1, 0.15, and 0.2 




60 80 

(GeV) 

Thermal Avg. av at x=24 by Final State, for mH=140 GeV, e=0.05 



SM fermion final states 
W*W" final state 
ZZ final state 
HZ final state 
HH final state 




40 60 80 100 

Mass of dark matter candidate (GeV) 



140 



Figure 3: Top: Relic density versus dark matter mass for Z suppression of (top to 
bottom curves) 0.05, 0.1, 0.15, and 0.2. Bottom: Thermally averaged annihilation 
cross section {av) by final state at x = 24: (near freeze-out) versus mass, with Z 
suppression of 0.05. Both plots assume a Higgs mass of mh = 140 GeV. 
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Relic Density for mH=160, e=0.05, 0.1, 0.15, and 0.2 
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Thermal Avg. av at x=24 by Final State, for mH=160 GeV, e=0.05 
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Figure 4: Top: Relic density versus dark matter mass for Z suppression of (top to 
bottom curves) 0.05, 0.1, 0.15, and 0.2. Bottom: Thermally averaged annihilation 
cross section {av) by final state at x = 24: (near freeze-out) versus mass, with Z 
suppression of 0.05. Both plots assume a Higgs mass of nih = 160 GeV. 
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based on current data. We wish to consider the radiative stabihty of such small e. 
The discrete symmetry under which 

ih,h^) ^ (/i^-Zi) (6) 
(Li,L2) ^ {-L2,L,) (7) 

(^1,^2) ^ (^2,^l) (8) 
(Sl,S2) (s2,Sl) (9) 

implies that the neutral mass eigenstates are maximally mixed Dirac neutrinos and 
that the charged states are degenerate. It is, of course, broken by Standard Model 
hypercharge and fermion mass differences. However, the Standard Model breaking 
does not generate a deviation from maximum mixing at one-loop order. We have not 
analyzed the two- loop contributions carefully, but estimate roughly that, assuming 
this symmetry is exactly preserved in the dark sector at Mq = 10^^ GeV, the Standard 
Model breaking generates an e ~ 0.001 — 0.01 at the electroweak scale. If the charged 
state Yukawa couplings are not exactly degenerate (but e = at Mq), then the 
non-degeneracy Yi = 1^2(1 + at Mq radiatively generates at the weak scale a Z 
suppression 

_ 3 ^JMa\Kl±jq. 



167r2 \mwj Kl-Kl ' 

where Ki and K2 are the eigenvalues of k. For reasonable values of the parameters, 
we find that e ~ 5. These estimates are justified in Appendix lA. 21 

One can also consider an approximate SU{2) symmetry that mixes the lepton 
doublets and rotates the Higgs. It acts as 

n = {-h',h)^ULm(R, (10) 

C = {LuL2)^UlCUr. (11) 

and can act on the Standard Model quark doublets as Q = {u, d) QUr. The 
symmetry is thus an extension of the ordinary Standard Model SU {2)r that enforces 
the custodial SU(2) p = 1 relation, though it acts very differently on the vector 
leptons than on the Standard Model fermions. The only invariant combination of Ti 
and C is 

rr[(iaf ^)7i:^(zfTf ^)/:] = hLi + h'L2 (12) 

Unlike the discrete symmetries above, it is broken by Yi and Y2 even when they are 
degenerate, as well as by K unless Ku = —K2i (i.e. 6^ = 7r/4, with x = 75(^1 + ^2) 
massless). In the limit that this symmetry is exact in the neutral sector, the mass 
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matrix takes the form 



V~ll/V O/V — sinQ;(//) cosq;(//) 

The condition = is equivalent to demanding that some hnear combination 
So = — sin asi + cos as2 of the singlets is non-interacting, and so a mass for x will 
never be generated radiatively within the SU{3) x SU{2)l x U{1)y effective theory. 
Small explicit breaking is needed and generically results in a Z-suppression e of the 
same order as m^/m^i where m^i is the mass of x' = ~ ^2)- 

4.2 An Adjoint/ Vector Lepton Model 

Another minimal extension of the Standard Model that can naturally implement the 

electroweak dark matter scenario consists of a vector pair of weak doublets Li, L2, an 

2^0 '-p+ 



SU{2)i adjoint fermion T = ^_ J , and a singlet s (this is the particle con- 

tent of Split-SUSY without a heavy gluino, and with no trace of the supersymmetric 
relations). As in the vector lepton model, we impose a chiral symmetry to suppress the 
invariant mass terms and mixing with the Standard Model. The Yukawa Lagrangian 
is 

dnt - Yi/iTLi + Y2N'TL2 + Y^hL^ + ¥^^"12. (14) 
An approximate SU {2)r symmetry acting as 

n = {-h',h) ^UlHUr (15) 

L = {L^,L2)^UlCUr, (16) 
with T and s both singlets under SU{2)ji requires the Lagrangian to be of the form 

Cir^t = y'Tr[(iaf))H^(iaf))T'^]+^r^[(^^f^)^^(^t^f^)'^]5 

= Y'{hTLi + h'TL2) + Ys{hLi + h''L2), (17) 

requiring Yi — Y2 — Y', Y^ — Y4 — Y. The corresponding low-energy mass matrix is 

= -(-''('0 :)(;») -S<-^—-^'' 
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where m = Yv, m' = Y'v, x = ^{Ni + N2), and x' = 71(^1 " ^2). Thus, SU{2)r 
generates maximally mixed neutrino mass eigenstates x x'- Standard Model 
hypercharge and isodoublet mass splitting break this symmetry, but communicate 
this breaking at two-loop. If all Yukawa terms are strongly coupled near the scale 
Mg ~ 10^^ GeV, then we obtain a low energy spectrum near ~ 100 — 140 GeV. 
Because the tree-level charged state masses are smaller than the x' mass, the dark 
matter must be the x ~ ^ state. To reproduce floM, we expect a x — s mass in the 
range of 35 — 45 GeV or 55 — 90 GeV. We have not studied the radiative effects in 
this model, but expect them to be of the same size as those discussed for the discrete 
symmetries of the previous section. 

4.3 Gauge Coupling Unification 

Both of the models discussed above can easily accommodate gauge unification, though 
this presents additional difficulties associated with proton decay or large contributions 
to precision electroweak observables. 

If we assume that there is not one vector lepton family but three (one for each 
Standard Model generation), then gauge couphng unification is achieved quite pre- 
cisely. Figure El shows two-loop gauge coupling unification predictions for as{Mz)- 
The predictions are within ^ 2a of the experimental uncertainty (the uncertainty in 
our calculation is ~ 0.004), and are in fact slightly closer to the observed a'^^^{Mz) = 
0.119±0.002 than in minimal SU(5) SUSY GUTs that predict a,(Mz) = 0.130±0.004 
(neglecting thresholds) [121 . However, the predicted unification scale is Mq ~ 10^^ 
GeV, so that avoiding dangerous proton decay rates is difficult in GUT completions 
of this model. This model also results in a large contribution to the precision elec- 
troweak observable S. With three families we have ~ ^, which is nearly 3-a from 
current best fits (see section 13^ . 

Adding a heavy SU(3) adjoint fermion (a "gluino") to the adjoint /vector lepton 
model recovers the low-energy spectrum of Spht-SUSY [T71 HH] , and unification is 
achieved quite precisely with a gluino mass of m ~ 8 TeV and a unification scale of 
Mg ~ 10^^ GeV. As discussed in Sec 15. 4| this also results in a large contribution to 
the S parameter. 

5 Experimental Prospects and Constraints 

In this section we discuss the implications of the dark matter scenario of the preceding 
sections for direct and indirect dark matter searches, collider searches, and precision 
electroweak observables. Figures IHl El and |H1 summarize experimental constraints on 
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Figure 5: Left: The top line (in red) shows the unification prediction for a^lMz) 
for the three family vector lepton model while the lower line (in green) shows the 
prediction for a two family model. The 1-a band for the experimentally measured 
as^Mz) = 0.119 ± 0.002 [2^ is shown by the dotted lines. Right: The top line (in 
green) shows the unification scale Mq for a two family vector lepton model while the 
bottom line (in red) shows Mq for a three family model. 



this model for run = 115, 140, and 160 GeV. We display the la and 2a regions con- 
sistent with nh"^ = 0.111 ± 0.006 (colored bands), the regions of parameter space 
excluded by direct Z-pole observations at LEP (upper left light green and blue curves 
are 3- and 2-a exclusions), and the constraints from direct dark matter searches (pro- 
cessed as described below). The dashed gray contour is the boundary of the DAMA 
95% CL. allowed region |^, while the three lower curves represent the EDELWEISS 
120! (dotted dark green), CRESST |2lj (dot-dashed black), and CDMS [22^ (dashed 
red) 90% excluded regions. Anything above the colored bands corresponds to a relic 
density below the observed total abundance, i.e. a heavy neutrino comprising only a 
fraction of the dark matter. 

The bounds from CDMS imply e < 0.01. There are three mass ranges for which 
the dark matter is not overproduced in this regime: at rriDM ~ 45, where even weakly 
coupled particles annihilate effectively through on-shell Z's, at 90-95 GeV, and near 
the Higgs resonance. In the latter two regions, the dark matter annihilates primarily 
through the Higgs channel. 

5.1 Direct Dark Matter Searches 

Direct searches for dark matter look for elastic and inelastic collisions between WIMPs 
in the local dark matter halo and heavy nuclei in detectors. Current searches are 
most sensitive to spin-independent interactions (which are generally mediated by 
scalar and vector nucleon currents) and probe the cross sections of these interactions 
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Constraints on DM mass and e (m|_|=1 15 GeV) 
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Figure 6: Exclusion plot in the e-moM plane for a Higgs mass of mn = 115 GeV. 
The purple and blue bands correspond to 1- and 2-a bounds on flnbm = 0.111 ±0.006 
from WMAP+CBI+ACBAR+2dFGRS data 0. The blue and green curves on the 
upper left correspond to the regions of partial Z widths of 0.2 and 1.7 MeV, cor- 
responding to 2- and 3-0" exclusions from the LEP Z-pole measurements, respec- 
tively j2]. The gray dotted contour is the 3-cr region of the DAMA signal ^^1; while 
the remaining curves are the most recent published 90% CL. excluded regions from 
EDELWEISS CRESST gl], and CDMS [2^ (from top to bottom, respectively). 
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Figure 7: The same as in Figure El for a Higgs mass of mn = 140 GeV 
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Figure 8: The same as in Figure IHl for a Higgs mass of mn = 160 GeV 
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using very heavy nuclei such as Ge, I, W or Xe 1211 121] • Current and future 
experiments probing spin-independent and spin- dependent WIMP-nucleon scattering 
include DAMA HSj, CDMS [221121, EDELWEISS 120], CRESST |2Z1I2II, PICASSO 
PI, SIMPLE [2S1, ZEPLIN jHU], XENON jSH, NAIAD [221, GENIUS ED and 
HDMS 

WIMP-nucleon scattering results are usually quoted assuming isospin invariance 
Gp ^ Gn, but for X the coupling to neutrons is Gn = while the proton coupling 
is greatly suppressed: Gp = — (1 — 4 sin^(^u;))G„ ~ — 0.076G„ [23j. We naively correct 
for this by reducing quoted sensitivities by ^(dz^)^°i^Z^ j = 0.28, 0.30, and 0.32 for 
Ge, I, and W detectors respectively [SEj- This reduction is taken into account in the 
exclusion plots El through |Hl but the cross-sections in the text are as quoted from the 
original papers, without the correction for isospin-dependent couplings. 

Although to date there has been no conclusive discovery of a WIMP signal, the 
DAMA experiment (a Na/I detector) has reported a 6.3a annually modulating signal 
attributable to a WIMP with mass m/jA/ ~ 50 GeV and WIMP-nucleon cross section 
cr ^ 7 X 10~^ pb [ini EZj. A mixed heavy neutrino with e ^ 0.02 — 0.04 would have 
a cross-section a^-n ~ 10^^ pb consistent with the DAMA signal. We see also that 
there are several masses {rriDM ~ 40 or 50 GeV) close to the reported DAMA signal 
mass for which this coupling would reproduce the observed ^Idm- 

The DAMA signal seems excluded at the 99.8% C.L. by the Ge detectors EDEL- 
WEISS and CDMS, which use stronger background discrimination techniques |j21 
I2ni 120] • The bound on WIMP-nucleon cross-section obtained from the CDMS null 
signal with standard estimates of the halo density and local halo profile is near 
cwiMP-n ^ 4 X 10~^ pb at m£)Af = 50 Gev (assuming isospin-invariance) j22]- This 
is well below the cross-section of cr^^n ~ 2 x 10^^ pb for a Standard Model heavy 
neutrino. Our model is consistent with this CDMS limit if e < 0.01 and with the 
EDELWEISS and CRESST limits if e < 0.02. If the one-loop radiative generation of 
e is dominant, we would expect 0.01 < e < 0.1 (see Sec. 14.11) . and <7x-n ~ 10^^ — 10~^ 
pb. Current experiments are probing this regime, and appear to favor smaller e and 
hence charged states that are very nearly degenerate. Two-loop effects from the Stan- 
dard Model generate 0.001 < e < 0.01 (10"^ < a^-n < few x 10"^ pb), and place a 
floor on the minimum e that can be obtained without resorting to fine-tuning. Cur- 
rent and future direct searches will increase sensitivities to awiMP-n ~ 10~^ pb, and 
so will detect a signal or rule out the viable parameter space for models of this type. 
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5.2 Indirect Searches 

Indirect signals of dark matter annihilations in the Earth, Sun or Galactic halo can 
also be used to probe dark matter. The decay products of annihilating WIMPs con- 
centrated by gravity can give rise to sizable fluxes of neutrinos, gamma rays and 
cosmic rays above the background of conventional astrophysical sources. We discuss 
here the general outlook for indirect searches in the electroweak dark matter sce- 
nario. We leave a detailed study of the indirect signatures of heavy maximally mixed 
neutrino dark matter for a future work. 

As a large body travels through the local dark matter halo, collisions between dark 
matter particles and nuclei in the body can dissipate enough energy that the dark 
matter particles become gravitationally bound and sink to the center of the body. An 
equilibrium is eventually expected between the rate of capture and of annihilations 
enhanced by the greater density of dark matter in the center of the body. Excess high- 
energy neutrino fluxes are the most reliable indicator of WIMP annihilations in the 
Sun or Earth, and result in bounds on spin-dependent and spin-independent WIMP- 
nucleon cross sections. Current and future neutrino telescopes carrying out such ob- 
servations include Super-Kamiokande [SHI, AMANDA [39 , BAIKAL [40j, ANTARES 
|41j . NESTOR |12] and IceCube To date, no statistically significant neutrino 

excesses from the Sun or Earth have been detected . 

In addition to neutrinos, gamma rays and high-energy cosmic rays from the galac- 
tic center also offer a strong indirect detection possibility jUj. Current and future 
experiments designed to measure gamma and cosmic ray fluxes include EGRET j45j . 
HEAT gB], PAMELA g7|, AMS gS], and GLAST 03. An excess in gamma rays 
above 1 GeV detected by EGRET HD] and in positrons peaking at 8 GeV by 
HEAT may be products of WIMP annihilations in the galactic halo jSUESESl- 
Moreover, a possible excess of microwave emission from the galactic center observed 
in WMAP microwave data is consistent with synchrotron radiation from e~^e~ pairs 
produced by WIMP annihilation 

The uncertainties in the annihilation rates and predicted flux excess due to uncer- 
tainties in the halo profile of the Milky Way make it difficult to reach any conclusion 
regarding the significance of these signals, but we consider briefiy their consistency 
with the model presented here. Previous authors have studied indirect signals com- 
ing from heavy stable 4th generation neutrinos [521 EHl EZl EH] ? but most recent work 
relating the observed gamma and cosmic ray excesses to WIMP annihilations have 
focused on the MSSM neutralino. The principal tension in explaining the HEAT 
and microwave signals from neutralino dark matter arises because, being Majorana 
fermions, their annihilations into light Dirac fermion final states such as e~^e~ are 
suppressed. As such, large bost factors are required to explain the strength of the 
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HEAT positron signal |S21 ESI- Moreover, as most of the electrons and positrons 
from neutralino annihilations are indirect products, neutralino models tend to pro- 
duce softer positron spectra, whereas the analysis of WMAP data in [S^ may prefer 
harder spectra. 

Our dark matter candidate is somewhat similar to the low-mass Kaluza-Klein 
state considered in |59j, in that both are Dirac particles. In that case, the HEAT 
signal could be explained with less need for large boost factors than in the neutralino 
scenario, and we expect the same to be true for our model. We further expect that, 
if mjjM ^ 80 GeV, our dark matter candidate would produce a harder spectrum 
of positrons than in the neutralino case, and perhaps be more consistent with the 
WMAP "haze" . Because our dark matter is Dirac, however, there may also be tension 
with the null neutrino flux signal from AMANDA |59J. Further work is necessary to 
compare the predicted gamma and positron spectra in our models to the EGRET, 
HEAT, and microwave observations. 

5.3 Collider Detection 

The most direct constraint on the dark matter particle itself, the lightest neutrino 
X, comes from the Z-width measurement at LEP. The measured partial width to 
invisible states is 499 ±1.5 MeV, obtained by subtracting the visible partial widths 
from the total width and assuming lepton universality. The standard model predicts 
a partial width of 167.29 ± 0.07 MeV to each species of neutrino, or 501.81 ± 0.13 
MeV overall [2]. If X is lighter than m^/2, then it will contribute a partial width. 



suppressed because of the reduced Z-coupling and a threshold suppression due to x's 
finite mass. Regions of exclusion at two- and three-cr {Tz{xx) = 0.2 and 1.7 MeV, 
respectively) are shown on the left side of the exclusion plots IHl El and |H1 It should 
be noted that the prediction of Tinv from the Standard Model alone is already 1.9o" 
above the measured width. 

A further constraint on this framework comes from the production of xx' pairs 
through an off-shell Z at LEPII. Though the XX^ and x'x'^ couplings are suppressed 
by e, the mixed couplings xx'^ are proportional to \/l — e^, and hence essentially 
unsuppressed. The maximum center of mass energy reached at LEP II was ~ 200 
GeV, so xx' pairs would be produced if + m^/ < 200 GeV. The x' would decay 
in the detector to a stable x and two jets or two lepton tracks via an intermediate Z. 
Thus, we believe that a two-jet or two-lepton track plus missing energy signal would 
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have been seen at LEP II if the above mass condition were met. We are not aware of 
constraints in the literature on this decay mode, and a more careful analysis of the 
statistics is required to determine the precise exclusion range. Because it depends on 

+ m^i, this bound is only constraining when combined with the assumption of 
high-scale perturbativity of the Yukawa couplings, which requires m^/ < 150 GeV. 

Because the dark sector interacts only weakly, and most of the energy in pp col- 
lisions at the LHC will be in the form of energetic gluons, production at the LHC 
will be accompanied by considerable backgrounds. The favored masses m^Af ^ 90 
GeV for the dark matter x and < 200 GeV for the heavier neutral and charged states 
make these models accessible to the LHC, but it is unclear how long it will take to 
obtain sufficient statistics to learn anything concrete about a given model pUj . 

5.4 Precision Electroweak Observables 

Any extension of the Standard Model that couples to the SU{2)l x U{1)y sector 
contributes to the gauge boson self-energies and hence modifies predictions for the 
oblique correction parameters S", T, and U [HI]. These contributions are particularly 
large for extensions involving fermions that get mass only from electroweak symmetry 
breaking, which do not decouple from the Standard Model even when irtf ^ mz- 

The T parameter measures the amount of custodial SU(2) breaking that occurs, 
(i.e. deviations from p = 1), and current constraints imply that Xli ~ (85)^ 

GeV^ at 95% CL, where Am^ = m? + — In — and the sum is over all 

isodoublets 0. SU(2)l fermion doublets also generically contribute positively to the 
5 parameter. Current fits suggest 5 = -0.13±0.10, T = -0.17±0.12, f/ = 0.22±0.13 

m 

The contribution of a single-family vector lepton model to 5* is 5* ~ which is 
sufficiently small to not be strongly at odds with precision electroweak data. More- 
over, under the assumption that the vector lepton sector is close to strong coupling 
near the GUT scale {Mq ~ 10^^ GeV in these models), the spectrum is not generi- 
cally split by enough to contribute largely to T. Only in cases where there is a light 
dark matter candidate (such as when the SU{2)^ symmetry of section l^?T] is approx- 
imately preserved) do we expect a sizable positive contribution to T. In this case, 
the experimentally preferred value of S is closer to 0, and hence more consistent with 
the positive contribution to S. Further analysis is needed to work out the detailed 
predictions for electroweak observables with a vector lepton family, but this minimal 
model is not ruled out. 

For the adjoint /vector lepton model, however, contributions to 5* are much larger, 
S* ~ -. As with the vector model, the contributions to T depend on the detailed 
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form of the spectrum but do not have to be large. Based on the S parameter alone, 
this minimal model is ruled out at ~ Aa (if there is a large contribution to T, then 
the exclusion is at ~ 3a). An invariant mass term for the adjoint field can 
be added while maintaining a dark matter mass arising entirely from electroweak 
symmetry breaking so long as my < Mz, thereby suppressing the adjoint field's large 
contribution to S. 

6 The Cosmological Constant, Gauge Hierarchy, 
and the Structure Principle 

In this section, we discuss some of the implications of the "structure principle" as 
defined in 0. For our purposes, this principle requires that large-scale structure 
develop in the universe as it cools. In the context of landscape scenarios, this principle 
is justified by the weak anthropic argument that biological creatures will not develop 
in a universe that cannot support the development of stars and other large-scale stable 
structure. We review how this principle was first applied by Weinberg to explain the 
smallness of the cosmological constant and discuss an apphcation suggested in [Hj to 
a possible resolution to the gauge hierarchy problem in the vector lepton model. 

6.1 Predicting the cosmological constant. 

Weinberg's argument predicting a cosmological constant that is very small if nonzero 
begins with the observation that gravitational collapse via Jeans instabilities can occur 
only after the universe becomes matter dominated. Moreover, linear sub-horizon 
perturbations to the energy density can only grow when the energy density in the form 
of a cosmological constant is smaller than the energy density of matter, so the universe 
must be in this regime when non-linear structures such as galaxies start to form. After 
matter-radiation equality, initial perturbations — scale as the acceleration parameter 
a, so non-linear structures begin to form after the universe has expanded by an 
amount (^)~^ after matter-radiation equality, and we require for structure formation 
A < Pmr{,—Y, where pmr is the energy density at matter-radiation equality. 

As was pointed out in jH], if the dark matter is dominated by cold relics a standard 
perturbative freeze-out calculation gives 




(21) 



where {av) is the thermal average of the annihilation rate and g^, is the number of 
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effective degrees of freedom. Weinberg's argument tlien bounds A* as, 

Ai/4 < 1^ til (22) 

Witliout any additional reason for A 3 to be small, this bound should be roughly 
saturated, so it is a generic prediction of the structure principle. 

For heavy (i.e. more massive than Mz) weakly interacting CDM particles, the 

2 

annihilation cross section av ~ ^ . The above argument then implies 



DM 



(Sp\3/4: 2 

^1/4 < (23) 

as pointed out in Empirically, heavy WIMP scenarios lead to ~ .l{ "^'^y )~'^ , 
thereby implying that m^M is of order a TeV and that A^/"^ ~ 'W^- 

Although weakly interacting TeV-scale dark matter seems to be empirically con- 
sistent, there is no a priori reason why rripiM is so much smaller than Mpi or so close 
to V. One simple possibility is that weak scale supersymmetry implies that ttidm ~ v 
as would also be the case for any other solution to the hierarchy problem that also 
contains a dark matter candidate. If we give up the assumption that the solution to 
the gauge hierarchy problem also explains dark matter, a simpler possibility emerges: 
// the dark matter particle gets mass only via electroweak symmetry breaking, then 
rriDM will naturally be close to v, and A^^'^ ~ is predicted! 

6.2 Generating the gauge hierarchy with the structure prin- 
ciple 

In addition to providing an explanation for the smallness of the cosmological constant, 
the structure principle also suggests a possible explanation of the gauge hierarchy in 
any model with a Standard Model Higgs and a dark matter candidate that gets mass 
only via electro-weak symmetry breaking, such as the vector lepton model discussed 
here. As we will show, electroweak symmetry breaking via a negative rn^ much below 
MpL now becomes essential for structure formation. If the Higgs mass parameter 
can be scanned in a landscape scenario and electroweak symmetry breaking can only 
happen in a very small window about its current value, then we will have found a 
possible explanation for the hierarchy. 

To see how this might work, we consider a suggestion first made in In the 
SM with a large top Yukawa coupling to the Higgs, the quartic coupling A has a 
negative UV fixed point Xuv ~ —yt and an IR fixed point Xip ~ +yf. If the Higgs 
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is sufficiently light (i.e. A is small), then for energies beyond some threshold Mcross, 
A will be driven negative by RG running and the vacuum becomes meta-stable. If 
|A(//)| does not become too large in the UV, the vacuum in the early universe has only 
a very small amplitude to decay as the universe cools, and the theory is safe from 
cosmo logical disaster. Turning this picture around, suppose A(M*) starts negative at 
the cutoff scale, say A(M=k) ^ Xuv- As the universe cools A(/i) crosses zero at a scale 
Mcross exponentially suppressed with respect to the cutoffi Moreover, the low energy 
universe looks very different depending on m'j^. 

Case I {m'jj > 0): The higgs vacuum centered around zero is exactly stable and 
electro-weak symmetry is broken only by the QCD quark condensate. Because M^ ?a 
«2Aqcd, sphaleron transitions operate all the way down to Aqcd, biasing zero baryon 
number and washing out any baryon asymmetry in the universe down to the freeze- 
out level of ^ ~ 10^^^. (See jU] for a discussion of this effect in more detail). 
The universe remains essentially radiation dominated and devoid of baryons, thereby 
prohibiting structure formation. Thus, the structure principle would rule this range 
of parameters out. 

Case II {mjj- < —M^^^^^): The Higgs field is heavier than Mcross, so Xijan) < and 
the Higgs vacuum is unstable. Fluctuations trigger (0) — > Mpi. When the vacuum 
decays, the universe becomes dominated by a cosmological constant A ~ Mpj^, so no 
structure can form in this case either. 

Case HI {—Mcross < < 0): Becausc electroweak symmetry is broken below 
the scale Mcross, aH of the matter fields can obtain sizable masses and the universe 
can become matter-dominated allowing structure to form. 

Thus, if the dark matter gets mass from electroweak symmetry breaking, the 
prediction without any fine-tunings is that the Higgs mass should be close to Mcross- 
To determine how close should be to Mcross to avoid fine-tuning, we calculate 
the 1-loop effective potential and look for the largest values of \m'jj\ for which a meta- 
stable minimum develops in l^((/>)i_;oop at nonzero (j). Following p], we approximate 
A(/i) near Mcross using an approximate solution to its RGE. Following appendix O 
we see that A(u) ~ —b\og(YT^ — ) where b is an RG coefficient and, 

YMlzl^ ^ -e^(^) - -log(^)(^)^ (24) 
^M ' 2 W '^M ' ^ ' 

cross ^^^cross ^ ^^^cross ^^^cross 

where e = \mH\/ Mcross- For the Standard Model, h ^ 0.076 and numerical solutions 
indicate that a stable secondary minimum develops for e < 0.09 at (0) < {]AMcross- 
Typical values for slightly smaller e are in the range of {(f)) < Q.2Mcross- Thus, without 
fine-tuning to a scale below that suggested by the structure principle, we would expect 
physical Higgs masses in the range of (0.1 — Q.2)Mcross for the Standard Model. With 
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Figure 9: Mcross versus heavy degenerate charged and neutral state masses. The one- 
loop running of A was computed using tree level values of the heavy lepton yukawa 
couplings. We used physical Higgs masses of mn = 155, 135 and 115 GeV from top 
to bottom and a top mass of = 178.1 GeV. 



a more careful analysis of the one-loop effective potential, we expect to find that the 
threshold mass is slightly higher than this. 

In the Standard Model with a higgs mass of 115 GeV, we calculated M^toss for 
a top mass in the experimentally allowed range of 169.2 — 188.5 GeV. Mcross ~ 
7, 40, 120, 1350 TeV for top masses of = 188.5, 178.1, 174.3, 169.2 GeV respectively. 
If the top were mt ~ 188 GeV, then a higgs mass of 115 GeV would be entirely 
consistent with this scenario without fine-tuning. A top mass near its central value 
of = 178 leaves a small hierarchy to deal with. 

Turning to the vector lepton model, the additional Yukawa couplings drive A(/i) 
negative even faster. Thus, we expect to find Mcross lower for heavier lepton masses. 
Setting all heavy lepton masses equal to My, we consider the dependence on My of 
Mcross for several values of Higgs mass in figure IHl The hierarchy between typical 
values of Mcross the range of 10-50 TeV in the Standard Model is eliminated for 
a light Higgs and a modestly heavy vector lepton spectrum of ~ 150 GeV, but A(/i) 
runs dangerously negative in the UV. 

In the Standard Model, requiring that A(/i) not run so negative that the vacuum 
should have decayed during the last 10^° years leads to a bound on the Higgs mass 
of mH > 115 GeV. For the vector model extension, even with a light spectrum near 
~ 100 GeV, A(/i) tends to run too negative in the UV. In our analysis, we calculated 
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Figure 10: The vacuum instability scale Mdecay versus heavy degenerate charged 
and neutral state masses. The one-loop running of A was computed using tree level 
values of the heavy lepton yukawa couplings. We used physical Higgs masses of 
mfj = 155, 135 and 115 GeV from top to bottom and a top mass of = 178.1 GeV. 



the scale at which A(/i) becomes sufficiently negative for the vacuum in our universe 
to have decayed already (assuming mn > 115 GeV). We discuss the details of this 
requirement in appendix O For this rough analysis, we require that A(yu) not be less 
than Xdecay ~ —0.13. Figure Uni displays our results. We found that requiring A to not 
run past the stability bound requires mfj > 155 GeV. Thus, for the above explanation 
of the gauge hierarchy to work, new physics must enter below the scale Mdecay to 
prevent A(/i) from running too negative. This new physics should make the UV fixed 
point for A less negative. For example, new SU{3) fermions at an intermediate scale 
slows the running of ^3 and hence helps keep the top Yukawa yt small, in turn making 
the UV fixed point for A less negative. Alternatively, gauging the chiral symmetry 
that was introduced to forbid mass terms with a strongly coupled U(l) that is broken 
near Mdecay (which requires adding an extra vector lepton generation to eliminate 
anomalies) drives the charged and neutral state Yukawa couplings smaller in the UV, 
increasing the UV fixed point of A. 

Finally, the upper bound derived from the requirement of not hitting any Landau 
pole up to Mg ~ 10^^ GeV requires mn < 165 GeV. The exact limit depends on 
vector lepton masses. Consequently, the self-consistency region for the Higgs mass in 
our minimal vector lepton model is 155 < ^ 165 GeV. 
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7 Conclusions 

We have considered a simple and well-motivated explanation for the origin of dark 
matter, namely that it consists of particles that get their mass entirely through elec- 
troweak symmetry breaking. This framework connects the scale and existence of 
massive dark matter quite simply to the electroweak scale. Moreover, this explicit 
connection affords new possibilities for explaining the smallness of the cosmological 
constant in relation to the weak scale as well as the gauge hierarchy. 

In the simplest models studied here, the dark matter candidate is a mixture of 
two Dirac neutrinos with opposite isospin, and so has a Z coupling suppressed by 
e = cos(26'). Several approximate discrete and continuous symmetries can guarantee 
the radiative stability of e 0. If e < 0.15 is generated radiatively at one-loop 
order from breaking of the approximate symmetry within the dark matter sector, the 
observed relic density can be accounted for with a dark matter mass in the range 
of 30 — 35 GeV or 55 — 95 GeV. This is, however, disfavored by the null results of 
several direct WIMP searches. The range e ~ 0.01 — 0.001 is typical if the approximate 
symmetry is broken only by the Standard Model and e is generated at two-loop order. 
Then the dark matter mass is predicted to be ttidm ~ 45 GeV or moM ~ 90 — 95 
GeV and WIMP- neutron spin-independent cross sections are a wim p-n ~ 10-^-10-^ 
pb. Current and future experiments sensitive to awiMP-n ~ 10~^ pb will probe this 
regime in the next few years and either detect a signal or rule out the viable parameter 
space of this class of models. 

We found that gauge couplings in a three-family vector lepton model unify quite 
well near Mq ~ 10^^ GeV. A second model with the particle content of Split-SUSY 
unifies precisely at Mg with a heavy "gluino" mass of ~ 8 TeV. Future work is needed 
to build realistic "top-down" GUT models. Further work is also needed to examine 
consistency of the gamma ray and positron spectra across the ~ 1 — 100 GeV range 
expected from WIMP annihilations in the galactic halo with HEAT and EGRET 
results, as well as to consider other indirect search possibihties. 
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A Renormalization Group Results For Vector Lep- 
ton Model 

A.l Vector Lepton Model RG Equations 

We have used the following /3-functions in Sections 4 and 6. They are derived from the 
general expressions in [021 ESI El]- Here t = In/i, where /i is the renormalization scale. 
gi = \J\g' is the coupling of the S'f/(5)-normalized U{1)y- We quote two-loop beta 
functions for the gauge couplings, one-loop for the Yukawa and quartic couplings. 

1. Gauge couplings (2-loop) 



(4vr) 




df{Y,^ + Yi)-dftriK^K) 



(25) 



with 
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2. Yukawa couplings (1-loop): 



, . o (Iitx/\^ f , 9o 9o 3o 3o o9 ^9 \ o 

(47r) = fi;i2 [y-iKS] - —g^ - -g^ + -k^^ + -^^2 + 3^22 - -Y^ J + 3fi;nfi;2ifi;22 

, , (-) G?Av21 iT7-/r^\ ^9 ^9 9 ^9 ^9 ^t^9\ 

(47r) ^ = K21 ( ^^l-^) - - -^92 + + -^21 + 2'^22 - ) + 3K12«12K22 

, \ 9 ( f (-^\ ^2 ^2 o2 ^2 ^2 7"2 \ o 

(47r) ^ = K21 ( >2('5) - —^1 - -^2 + 3^12 + 2'«2i + 2'^22 - 0^2 ) + 3«;i2Kl2«21 

^dYx _ ^ ( I n\ 9 „2 9 2 3 2 3 2,3 2 



(4^f^ = Y,{Y,(S)-\gl-\gl-\^l,-\^,, + \Y:i^, (26) 
with 

y2(-S) = Tr (3y„ty^ + gy^ty^ + y^ty^) + ^2 ^ y^2 ^ tr(ii:^K) 

where Tr is over the three Standard Model generations, and tr over the two 
dimensions of K. 

3. Higgs quartic (1-loop) 

+ AY2{S)X-AH{S), (27) 

with 

H{S) = Tr (3(y„ty^)2 + 3(y^ty^)2 + (y^ty^)2) ^ ^4 + ^^4 ^ tr((i^^i^)2) 

A. 2 Radiative Generation Of e 

We note first that for small A, the matrix 



fi;/y2(l + A) /tVV2(l + A) 
-«;/^(l-A) fi;/^/2(l-A) 



(28) 



can be diagonahzed by a rotation of ^ = 7r/4 + ^^^^A on the left and ^^§^^A 
on the right. The maximally mixed Yukawa matrix is of this form with A = 0, 
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and all radiative breaking effects we will consider can be parameterized as radiative 
generation of nonzero A. Thus, we can immediately relate a splitting A of the Yukawa 
matrix elements to the associated Z-coupling suppression e = cos(2^) = 2^^3^A. 

For the estimates that follow, we work in the approximation that, although the 
eigenvalues k and k' of K also run, their fractional change is small enough that it 
does not greatly alter the running of A. 

We consider first the radiative effect of a splitting Yi = ¥2(1 + S) of the charged 
state Yukawa couplings at Mq, assuming maximally mixed K at Mq- From the RGEs 
above, 

so an order-of-magnitude estimate of e at the weak scale is given by 

eiMw) ~ Hjf )T7-2-^^Y2S ~ 6 (30) 

assuming k'/k, ^ 1.5 — 2. 

In the case where Yi = I2 at Mg is exact, e is not generated radiatively at one- 
loop, but there is a two-loop contribution to /9(A) from the Standard Model that goes 
as yfg2- Naively, we expect 

esAiiMz) - 7-^ ^2 ^ '^,2 ^ta2N, (31) 

where N accounts for color factors, the prefactor of the diagram in the RGEs, and the 
factor of 2 in the running of e. We estimate a radiatively generated e ~ 0.001, though 
this could easily change by an order of magnitude if is large or small, or because 
of the running of the parameters. Because this radiative generation appears to put 
e in the region of parameter space preferred to avoid direct detection, this radiative 
generation merits more careful study. 



B Calculating Relic Dark Matter Abundances 

In this appendix, we review our calculation of the relic density for the heavy stable 
state X of the vector lepton model of section 14.11 Our calculation is based on the 
discussion in [12], [Mj, and 

Consider the evolution of the number density ni and n2 for x (with mass m)and 
its antiparticle x- Making the standard assumptions of I3j and neglecting the co- 
annihilation effects discussed in ^3]; the evolution of the total number density n = 
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ni + n2 is described by the Boltzmann equation, 

n + 3Hn = ^{avM){n^ -n%), (32) 

where Ueq is the equihbrium number density, the Moller velocity vm is defined so that 
fM^i^2 is Lorentz invariant, and {ctvm) is the thermal average of the annihilation 
cross section calculated using the numerical integral 

1 f°° 

= Sm^TKiim/T) J,^. ' 4-^)v^^i(v^/^)^^' (33) 

where Ki are modified Bessel functions of order i. It is useful to treat the effect of the 
expansion of the universe implicitly by using the variable Y = ^ where s is the total 
entropy density of the universe. We also use a dimensionless temperature variable 



a; = ^, so that IHzl becomes 

Y = -s{avM)iY'-Y,l). (34) 

From this we can obtain, 

dY 1 ds 



dx 3H dx 



{avM){Y'-Y'), (35) 



where the Hubble parameter is H = y^nGp with G the gravitational constant and 

p the energy density of the universe. Using conventional definitions p = geff(T)^T^, 
and s = heff(T)^T^ jTS], we can finally re-cast the Boltzmann equation into the 
form 



where 



% = H^r^''-^{-VM)iY^ - ny, (36) 

gj/f ^heff dT 

In our analysis, (ctwm) was calculated numerically using equation ESI and then 
used to numerically solve the Boltzmann equation 1221 Throughout, we assumed that 
the heavy charged states all had mass ruheavy = 140 GeV and the heavy neutral state 
m^/ ^ 1.3 ~ consistent with a GUT-scale strong coupling assumption and sufficient 
for coannihilations to be negligible |14U15j . The annihilation channels included in our 
analysis are listed in Tabled Feyncalc was used to simplify the trace and Lorentz 
algebra involved in calculating the Feynman diagrams. Mathematica version 5.1 was 
used for all numerical calculations. 
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Final State 


s- channel 


t-channel 


u-channel 


ff 


Z, h 








Z Z 


h 


X) X 


X) X 


w+w- 


h, Z 


El 




Z h 


Z 


X 


X 


h h 


h 


X 


X 



Table 1: List of XX Final State annihilation channels included in our relic den- 
sity analysis. The entries in the table are the intermediate states contributing to 
annihilation into a given final state for a given channel. 



C Higgs Vacuum Stability 

In this appendix we briefly review our calculation of the Higgs vacuum stability bound 
in the vector lepton model of section 14.11 The physics of vacuum decay in relativistic 
quantum field theories is discussed in jSHl IHZ] and an analysis applied to the Standard 
Model can be found in 

With the addition of a single vector lepton family coupled to the Standard Model 
as in section ing the Higgs quartic running is modified to include yukawa contributions 
from the heavy charged and neutral states, 

{infPx = 12A2 + 4(3y2 + 2yl + - 4{3yt + 2y% + «:^). (38) 

For a small Higgs mass and hence small X{mh), the large Yukawa couplings yt and ye 
will drive A(yu) negative for large n, thereby making the vacuum unstable. We wish to 
demand that in the lifetime of the universe and neglecting finite-temperature effects, 
the probability for the vacuum to decay is less than 1. This requirement furnishes a 
lower bound on A. Applying the results of [UH], the decay rate per unit volume for 
the nucleation of a bubble of true vacuum of radius i? is ~ (-^)^e~'^^, where Se is 
the Euclidean action of the "bounce" solution (we are dropping multiplicative factors 
of order unity). Using Se = 3X^7^) we take the probability for vacuum decay pd to 
be, 

1 f dRflV 16.^ , , 

771/ If (5 <^9) 

where H = 2.133/i x IQ-^^ ^ 15x4 ^ 10-^2 GeV {h = 0.71) is the Hubble constant 
and accounts for the spatial and temporal extent of the universe. For the class of 
models under study, A runs negative and then tracks its approximate UV fixed point. 
Thus, in practice the integral ()39|) is dominated at the scale where A is most negative, 
Mnuc (this is typically the GUT scale Mq ^ 10" - 10^^ GeV in these models). In 
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this case, we take, pd ~ MqB ^^<-'^g) to get 

|A(M„)| < (40) 
Using Mg ^ 10^^ GeV, we find Xdecay ~ -0.13. 
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